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Abstract 

Due to the AdS/CFT correspondence the question of instability of Anti-de-Sitter spacetimes 
sits in the intersection of mathematical and numerical relativity, string theory, field theory and 
condensed matter physics. In this essay we revisit that important question emphasizing the power 
of spectral methods and highlighting the effectiveness of standard techniques for studying nonlinear 
dynamical systems. In particular we display explicitly how the problem can be modeled as a 
system on nonlinearly coupled harmonic oscillators. We highlight some of the many open questions 
that stem from this result and point out that a full understanding will necessarily required the 
interdisciplinary cooperation of various communities. 
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The Anti-de-Sitter/Conformal Field Theory correspondence (AdS/CFT) has been responsible for 
the growing interest in spacetimes with negative cosmological constant A < 0, that is, in asymptotically 
locally AdS spacetimes. The main content of the AdS/CFT correspondence is the identification of 
asymptotic data of fields in AdS with sources and vacuum expectation values of operators in the dual 
field theory M. 

The AdS/CFT correspondence was readily expanded to include field theories at finite tempera- 
ture. In this context a field theory in equilibrium at finite temperature is dual to a black hole in 
asymptotically AdS spacetime. One very important development has been the establishment of the 
correspondence at the level of small fluctuations. More precisely, fluctuations of systems at equilibrium 
or weakly driven near equilibrium satisfy a universal relation known as the fluctuation-dissipation the- 
orem which connects spontaneous fluctuations to the linear response. The search for similar relations 
for systems far from equilibrium has been an important problem for many years with applications 
in ion collisions at Relativistic Heavy Ion Collider (RHIC) and Large Hadron Collider (LHC) and 
also in condensed matter physics. The AdS/CFT correspondence has successfully incorporated those 
near-equilibrium relations and is now posed to take the next step and tackle far-from-equilibrium 
phenomena. To study far-from-equilibrium phenomena one needs to study the evolution of Einsteins 
equation with appropriate boundary conditions. Therefore the study of gravitational collapse has 
emerged as one of the most interesting problems in the interface of the interests various communi- 
ties including, among others, the numerical relativity community, string theory, high energy particle 
physics and condensed matter. 

An important property of asymptotically AdS spacetimes, which is different from our intuition in 
asymptotically flat spacetimes, is the presence of a timelike boundary at spatial and null infinity where 
suitable boundary conditions need to be prescribed. 

In this context the study of the dynamics of a self interacting massless scalar field in AdS spacetime 
has recently been initiated (see for example [Q, ^]). Even at this very incipient state very interesting 
results that remind us of the richness of the nonlinearities of Einstein's field equations have already 
been obtained. Before entering into the qualitative description of these features, let us consider the 
model studied by Bizon and Rostworowski (BR) § i n which the spacetime is described by, 



where I 2 = — 3/A, and the metric functions A, 8 depend on (t,x). Here t is the time coordinate and 
< x < tt/2 covers the spatial domain with the boundaries at the origin (x = 0) and the spatial 
infinity (x = 7r/2). To write down the field equations we define the auxiliary variables n and & as, 



where (p{t, x) is the scalar field, prime and dot means derivatives with respect to x and t, respectively. 
Now, the Einstein-Klein-Gordon system reads, 
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which show the evolution and constraint constraint. The most simple exact solution of this system 
is the pure AdS spacetime characterized byII = <I > = 0,5 = and A = 1. As a matter of fact this 
solution can be viewed as analogous to the Minkowski spacetime when A = 0. 

According to BR the scalar field evolves by performing a sequence of reflections between the 
boundaries and eventually is driven to instability to finally collapse forming a black hole. They have 
shown that this feature is valid for any arbitrarily small initial amplitude of the scalar field implying 
instability of AdS spacetime, contrary to the stability of the Minkowski spacetime with respect to 
arbitrarily small perturbations which can essentially dissipate to spatial infinity. The main reason 
behind such surprising result is the fact that the timelike nature of the boundaries of AdS spacetime 
act as a confining box. The continued reflections of the scalar field eventually excite modes with 
higher and higher energies that produces an inevitable road toward the instability suggesting typical 
signature of a turbulent system. Further supporting evidence to the claim of BR was provided in 
Dias, Horowitz and Santos § who considered the purely gravitational problem of nonlinear stability of 
AdS by studying non-spherically symmetric configurations. This work was performed perturbatively 
and therefore could not explore the full nonlinear regime such as horizon formation. Nevertheless it 
gathered evidence to the claim that energy is transferred to ever higher and higher frequency modes. 

In this essay we demonstrate the efficacy of spectral methods || in gaining understanding of the 
underlying physics. For this end we have solved numerically the field equations using spectral methods 
for which a brief description follows. The start up and also the central idea of any spectral method 
is to approximate the relevant fields - n, A, 5 - as appropriate series with respect to certain basis 
functions. For instance, taking the field IT, we have, 

N 

u a {t, y ) = J2 a k(t)My), (0-6) 

fc=0 

where dk(t) are the unknown modes, the variable y is a rescaling of x, y = 4x/tt — 1, so that the origin 
x = and the spatial infinity x = tt/2 are located at y = — 1,+1, respectively. The basis functions 
ipk(y) are constructed such that the conditions for smoothness at the origin, at spatial infinity and 
also the finiteness of the total mass-energy content of the spacetime are satisfied. This means that, 



U(t,y) = n (t) + O((y + l) 2 ), neary 
U(t,y) = 0((y-l)- 2 ), neary = 1. 



(0.7) 
(0.8) 



One of the main consequences of the spectral method is to transform a partial differential equation 
into a finite set of ordinary differential equations. In the present case these equations are for the modes 
modes dk(t). 

Following BR we consider initial data of the form, 



$(0,y) = 0, n(0,y) =e exp 



tan 2 (vr/4(l + y)) 



(0.9) 



where we have fixed a = 1/4 and eo is the initial amplitude of the scalar field distribution. As a 
very useful illustration of the reflections of the scalar field, we show in Fig. 1 a sequence of plots of 
Il(i, y) evaluated in distinct intervals of time. The first plot correspond to the initial distribution in 
which we have selected eo = 0.5 and the other plots are evaluated in several subsequent instants. The 
arrow indicates the sense of propagation that is toward the boundary y = 1 and after a reflection 
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returns to the boundary y = — 1. This overall periodic motion remains for a long time before the 
instabilities triggered by the gravitational sector begin to dominate the dynamics. The greater the 
initial amplitude the less time the system undergoes this periodic motion. 
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Figure 1: Sequence of snapshots showing the reflection of the scalar field starting from the initial data 
(first graph). The arrows indicate the sense of propagation and this sequence is about to repeat for a 
long time. Here the boundaries y = —1,1 represent the origin and the spatial infinity. Here eo = 0.5. 



Let us discuss a convenient way of viewing the dynamics of scalar field in AdS spacetimes. For 
the sake of simplicity in the weak field regime, 5, 1 — A <C 1, the Einstein-Klein-Gordon equation in 
spectral space can be written as, 

~~ I — 2i 7 w [tan 2 (x(y))]V.(t) +T(S,1 - A, a k , a k , a k , ■■) = 0. (0.10) 
tan [x{y)) 

In this way, the first two terms represent the evolution of the linearized modes (see Wald Q) of the 
scalar field that oscillates in well defined frequencies oj? = (3 -+- 2j) 2 , j = 0, 1, 2, and the remaining 
denoted by J- contains the interaction with the gravitational sector. Therefore, we imagine our system 
as a set of nonlinear interacting oscillators, and consequently the periodic motion can be understood 
as taken place on a torus in the abstract phase space spanned by the modes a k (t). Further numerical 
experiments indicated that if the term T is initially small, its influence on the long term dynamics 
is to alter the natural frequencies of oscillation producing a quasi-periodic motion. Then, the motion 
becomes unstable and the torus is eventually destroyed forming a strange attractor, which is a signature 
of chaos in our dynamical system. This means that the evolution of the scalar field in AdS is turbulent. 
This route to turbulence resembles the Ruelle-Takens scenario |7j where the initial amplitude eo plays 
the role of the control parameter. In Fig. 2, a visualization of this route by a sequence of plots of n 
and <p evaluated at the origin for three values of the initial amplitude, eo = 0.5, 3.0, 5.0. These curves 
can be understood as the projections of a trajectory in phase space into the plane (n, ip). In the first 
plot the closed curve means that the motion is periodic as expected since the initial scalar field is very 
weak (eo = 0.5). By increasing the amplitude the curve is no longer a closed one as a consequence of 
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Figure 2: The increase of the initial amplitude eo produces the transition from the regular motion on a 
torus (eo = 0.5) represented by a closed curve on the plane (H,(fi) to quasiperiodic motion (e=3.0), and 
to turbulent motion (eo = 5.0). In all cases tt = 30.0. 



the new frequencies producing a quasiperiodic motion, and finally for eo = 5.0 the structure is typical 
of a chaotic system in which the trajectory seems appear to skip around randomly. 

The last numerical experiment consists in determining the power spectrum of the scalar of curvature 
R = l// 2 [^cos(x) 2 (-n 2 + $ 2 ) - 12] evaluated at the origin. We basically have set eo = 0.001 in order 
to reproduce for a long time the linear behavior, and eo = 5.0. In the region of frequencies under 
consideration, it is clear the noisy structure of frequencies. Also, a broad band of high frequencies are 
also present typical of the turbulent regime. 

Many important questions remained unanswered and this only witnesses the fact that the dynam- 
ics of of Einstein's equations in asymptotically AdS spacetimes are substantially different from our 
intuition built on asymptotically Minkowski spacetime. 

The result that for an arbitrarily small perturbation the spacetime collapses after a few bounces 
challenges our intuition and leads to new venues of research into the understanding of the precise 
mechanism for collapse. 

Let us briefly discuss the some of the implications of this result for the dual large N strongly 
coupled field theories. There are many arguments leading to the fact that the gravity calculations 
presented here pertain to the microcanonical ensemble on the field theory side. First, conservation 
of energy tells us that we are in a system with a fixed amount of total energy. Further, since black 
holes in flat space have negative specific heat they can cool themselves by transferring energy to the 
surrounding and consequently be able to transfer more energy. We expect that small black holes in 
AdS will behave similarly. From the field theory point of view sometimes the canonical ensemble 
is relevant in which case it is worth revisiting the role of the Hawking-Page phase transition in the 
context of the AdS/CFT correspondence. The interpretation of the Hawking-Page transition in the 
context of the gauge gravity correspondence leads to the result that below a certain temperature, that 
is, below a certain radius of the black hole, the background that is relevant in the thermodynamic 
competition is no longer the Schwarzschild black hole but rather thermal AdS |§. This result means 
that for low temperatures the field theory goes into a confined phase rather than remaining in the 
quark gluon plasma. These suggests that the very low temperature black holes are not relevant from 
the field theory point of view as they correspond to an unstable phase when viewed in the canonical 
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Figure 3: Power spectrum of the scalar of curvature evaluated at the origin for two cases: for very small 
amplitude eq = 0.001, and for eo = 5.0. In the first case we have the regular or non-chaotic motion with 
frequencies expected from analytical considerations. In the second case the power spectral is typical of 
a turbulent system. The broadband frequency is present but not shown here. 



ensemble. This situations leads to the all important but conceptually difficult question of whether 
numerical relativity can accommodate a scenario that would correspond, on the field theory side, to 
the canonical ensemble. 

In this essay we have strived to show, by considering a particular example, that the study of gravity 
in asymptotically AdS spacetimes has many surprises and challenges in store for various communities. 
We have also tried to argue that the demand is also on methods. In particular, spectral methods seem 
more suited to answer questions about evolution of certain modes while finite difference methods will 
be more efficient in problems where a singularity forms. We have explicitly displayed techniques to 
study nonlinear dynamical systems with the end of understanding precisely the mechanism by which 
perturbations in AdS grow to the point where black hole formation occurs. 
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